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Abstract 

We investigate the phase diagram of the so-called Polyakov-Nambu- Jona-Lasinio model at finite 
temperature and non-zero chemical potential with three quark flavors. Chiral and deconfinement 
phase transitions are discussed and the relevant order- like parameters are analyzed. The results 
j_j ■ are compared with simple thermodynamic expectations and lattice data. We present the phase 

diagram in the (T, yUg) plane, paying special attention to the critical end point: as the strength 
of the flavor-mixing interaction becomes weaker, the critical end point moves to low temperatures 
and can even disappear. 

PACS numbers: ll.10.Wx, 11.30.Rd, 12.40.-y 



I. INTRODUCTION 

Symmetries play a fundamental role in physics. In fact, the modern fundamental physics 
is dominated by considerations underlying symmetries which can be exact, approximate 
(explicitly broken) and (spontaneously) broken. A special role is played by gauge or local 
symmetries which lead to the description of the "real world" as the so-called local/gauge 
theories with spontaneously broken symmetries. 

The concept of spontaneous symmetry breaking has been transferred from condensed 
matter physics to quantum field theory by Nambu [l|. It has been introduced in particle 
physics on the grounds of an analogy with the breaking of (electromagnetic) gauge symmetry 
in the theory of superconductivity by Bardeen, Cooper and Schrieffer (the so-called BCS 
theory). The application of spontaneous symmetry breaking to particle physics in the 1960s 
and successive years led to profound physical consequences and played a fundamental role 
in the edification of several models of elementary particles. 

Spontaneous breaking of chiral symmetry, in particular, is known to govern the low-energy 
properties of hadrons [2H4|. Some QCD-like models have been proposed before the advent 
of QCD, and the phenomenon of spontaneous breaking of chiral symmetry and Nambu- 
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Goldstone theorem was established more than 40 years ago [l]. 

The Nambu- Jona-Lasinio (NJL) model was proposed in 1961 to explain the origin of 
the nucleon mass with the help of spontaneous breaking of chiral symmetry 5|, |6[. At that 
time, the model was formulated in terms of nucleons, pions and scalar sigma mesons. The 
introduction of the quark degrees of freedom and the description of hadrons by Eguchi and 



Kikkawa 



BE 



in the chiral limit, where the bare quark mass is m = 0, and a more realistic 



version with m ^ by Volkov and Ebert 9Ml lj . initiate a very intensive activity by several 



research groups [12]. 



There is strong evidence that quantum chromodynamics (QCD) is the fundamental theory 
of strong interactions. Its basic constituents are quarks and gluons that are confined in 
hadronic matter. It is believed that at high temperatures or densities the hadronic matter 
should undergo a phase transition into a new state of matter, the quark-gluon plasma (QGP). 
A challenge of theoretical studies based on QCD is to predict the equation of state, the 
critical point and the nature of the phase transition. 

As the evolution of QCD at finite density/temperature is very complicated, QCD-like 



models, as for instance NJL type models, have been developed providing guidance and 
information relevant to observable experimental signs of deconfinement and QGP features. 
In fact, there has been great progress in the understanding of the properties of matter under 
extreme conditions of density and/or temperature, where the restoration of symmetries (e.g., 
the chiral symmetry) and the phenomenon of deconfinement should occur. These extreme 
conditions might be achieved in ultrarelativistic heavy-ion collisions or in the interior of 
neutron stars. In this context, increasing attention has been devoted to the study of the 
modification of particles propagating in a hot or dense medium [13|, |l4j. The possible 



survival of bound states in the deconfined phase of QCD 15H22J has also opened interesting 
scenarios for the identification of the relevant degrees of freedom in the vicinity of the phase 
transition 23N25|. Besides lattice calculations 26N29|. high temperature properties of QCD 
can be studied, starting from t 
the dimensional reduction 
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re^QCD Lagrangian, within different theoretical schemes, like 

34J. Actually 



3l| or the hard thermal loop approximation 
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both the above approaches rely on a separation of momentum scales which, strictly speaking, 
holds only in the weak coupling regime g ^ 1. Hence they cannot tell us anything about 
what happens in the vicinity of the phase transition. On the other hand, a system close to a 
phase transition is characterized by large correlation lengths (infinite in the case of a second 
order phase transition). Its behavior is mainly driven by the symmetries of the Lagrangian, 
rather than by the details of the microscopic interactions. 

Confinement and chiral symmetry breaking are two of the most important features of 
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36J have been successful 



QCD. As already referred, chiral models like the NJL model p,|6|, 
in explaining the dynamics of spontaneous breaking of chiral symmetry and its restoration 
at high temperatures and densities/chemical potentials. Recently, this and other type of 
models, together with an intense experimental activity, are underway to construct the phase 
diagram of QCD. 

The actual NJL type models describe interactions between constituent quarks, giving the 
correct chiral properties, and offers a simple and practical illustration of the basic mecha- 
nisms that drive the spontaneous breaking of chiral symmetry, a key feature of QCD in its 
low temperature and density phase. In order to take into account features of both chiral 



symmetry breaking and deconfinement, static degrees of freedom are introduced in t 
grangian through an effective gluon potential in terms of the Polyakov loop 



lis La- 
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44j. The 



coupling of the quarks to the Polyakov loop leads to the reduction of the weight of the 



quark degrees of freedom at low temperature as a consequence of the restoration of the Zjv c 
symmetry associated with the color confinement. 

In first approximation, the behavior of a system ruled by QCD is governed by the sym- 
metry properties of the Lagrangian, namely the (approximate) global symmetry SUz,(iV/) x 
SU n(Nf), which is spontaneously broken to SUy(iVj), and the (exact) SU c (iV c ) local color 
symmetry. Indeed, in NJL type models the mass of a constituent quark is directly related 
to the chiral condensate, which is the order parameter of the chiral phase transition and, 
hence, is non-vanishing at zero temperature and density. Here the system lives in the phase 
of spontaneously broken chiral symmetry: the strong interaction, by polarizing the vacuum 
and turning it into a condensate of quark- ant iquark pairs, transforms an initially point-like 
quark with its small bare mass rriQ into a massive quasiparticle with a finite size. Despite 
their widespread use, NJL models suffer a major shortcoming: the reduction to global (rather 
than local) color symmetry prevents quark confinement. 

On the other hand, in a non-abelian pure gauge theory, the Polyakov loop serves as an 
order parameter for the transition from the low temperature, TL^ a symmetric, confined phase 
(the active degrees of freedom being color-singlet states, the glueballs), to the high tempera- 
ture, deconfined phase (the active degrees of freedom being colored gluons), characterized by 
the spontaneous breaking of the Zat c (center of SU c (iV c )) symmetry. With the introduction 
of dynamical quarks, this symmetry breaking pattern is no longer exact: nevertheless it is 
still possible to distinguish a hadronic (confined) phase from a QGP (deconfined) one. 

In the PNJL model, quarks are coupled simultaneously to the chiral condensate and to 
the Polyakov loop: the model includes features of both chiral and Ztv c symmetry break- 
ing. The model has proven to be successful in reproducing lattice data concerning QCD 
thermodynamics [43|. The coupling to the Polyakov loop, resulting in a suppression of the 
unwanted quark contributions to the thermodynamics below the critical temperature, plays 
a fundamental role for the analysis of the critical behavior. 

One of the important features of the QCD phase diagram is the existence of a phase 
boundary in the (T, \xb) that separates the chirally broken hadronic phase from the chirally 
symmetric QGP phase. Arguments based on effective model calculations suggest that the 
QCD phase diagram can exhibit a tricritical point (TCP) /critical end point (CEP) where the 



line of first order matches that of second order/analytical crossover 45r|50||. The discussion 



about the existence and location of such critical points of QCD is a very important topic 



nowadays [51]. 



This paper is organized as follows. In Section [TT] we analyze the main features and 
symmetries of QCD. In Section HTT1 we present the model and formalism starting with the 
deduction of the self consistent equations. We also obtain the equations of state and the 
response functions. The regularization procedure used in the model calculations is also 
included. Section [IV] is devoted to the study of the equation of state at finite temperature. 
In Section [V] we study the phase diagram and the location of the critical end point. In 
Section IVIl we discuss the important role of the choice of the model parameters for the correct 
description of isentropic trajectories. In Section IVHI we analyze the effects of strangeness 
and anomaly strength on the location of the critical end point. In Section lYUH we proceed 
to study the size of the critical region around the critical end point and its consequences for 
the susceptibilities. Finally, concluding remarks are presented in Section IIXI 

II. SYMMETRIES OF QCD 

A. Quantum Chromodynamics 



The Lagrangian of QCD is written as 52|, [53] 



1 
4' 



£qcd = q{itD,-m)q - -jF^F? (1) 



where q is the quark field with six flavors (u, d, s, c, b, t), three colors (iV c = 3) and rh being 
the corresponding current quark mass matrix in flavor space (rh = diagy(m u , rrid, •••))• The 
covariant derivative, 

D, = 8, - igtAl (2) 

incorporates the color gauge field A a (a — 1, 2, ..., 8) and 

F% = %K~ duAl + gr hc A\Al (3) 

is the gluon field strength tensor; t a is the Gell-Mann color matrix in SU(3) 
([t a ,t fe ] = if a b c t c ,tr(t a t b ) = -^-) and f abc are the corresponding antisymmetric structure 
constants. Finally, g is the QCD coupling constant. 

The QCD Lagrangian is by construction symmetric under SU(3) gauge transformations 
in color space and because of the non-Abelian character of the gauge group QCD has some 



main features: it is a renormalizable quantum field theory [54| with a single coupling constant 
for both, the quark-gluon interactions and the gluonic self-couplings involving vertices with 
three and four gluons; it has confinement, i.e., objects carrying color like quarks and gluons 



do not exist as physical degrees o 
that has asymptotic freedom 55 



freedom in the vacuum. In addition, QCD is a theory 



561 ]. i.e., for large momenta, Q, or wavelengths of the 
order of 10" 1 fm (ultraviolet region) the couplings are weak and the quarks and gluons 
propagate almost freely. For low momenta, or wavelengths of about 1 fm (infrared region), 
it happens the opposite situation and the couplings are quite strong. The system is now 
highly non perturbative. According to this property, the attraction between two quarks 
grows indefinitely as they move away from each other. This implies that the interaction 
between quarks and gluons can not be treated perturbatively, making the perturbative 
treatment of QCD not applicable to describe hadrons with masses below ~2 GeV. 

The low energy regime is specially interesting once it is relevant to the study of hadronic 
properties as for example in low energy QCD and nuclear physics. 

The non-perturbative structure of the vacuum is characterized by the existence of quark 
condensates, i.e., it is expected non-zero values for the scalar density (qq), by the appearance 



of light pseudoscalar particles, who are identified with (quasi) Goldstone bosons [571, |58j , and 



also by the existence of gluon pairs [59] 



On the other hand, if QCD has a mechanism that includes the confinement of quarks 
the mass parameters m^ are not observable quantities. However, they can be estimated in 
terms of the masses of some hadronic observables through current algebra's methods. These 
masses are denominated by current quark masses to distinguish them from the constituent 
quark masses, which are effective masses generated by the spontaneous breaking of chiral 
symmetry in phenomenological models of quarks. 



B. Chiral Symmetry Breaking 

Due to the relevant role that spontaneous breaking of chiral symmetry plays in hadronic 
physics at low energies, this symmetry is one of most important symmetries of QCD. Here 
we will concentrate in the Nf = 3 case. In the chiral limit, i.e. m u = ma = m s = 0, QCD 
is chiral invariant which means that the QCD Lagrangian ([1]) is invariant under the group 
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Symmetry 


Transformation 


Current 


Name 


Manifestation in nature 


SLV(3) 

U V (1) 

SU A (3) 

Ua(1) 


q — > exp(-i^Hp-)g 

q — > ex.p(—iav)q 

q ->exp(-i 75A 2 ae ")q 

q -> exp(-i<y 5 a A )q 


V* = qi^q 

Vn = qi^q 

A% = qi^i^q 

A fM = qi^q 


Isospin 

Baryonic 

Quiral 

Axial 


Approximately conserved 

Conserved 

Spontaneously broken 

" Ua(1) problem" 



TABLE I. QCD symmetries in the chiral limit. 



of symmetries: 

U(3) L ® U(3)^=SUy(3) ® SU A (3) <g> U v (l) ® U A (1). (4) 

These symmetries are presented in Table [J where it is also possible to see the transfor- 
mations under which the Lagrangian is invariant, the currents that are conserved according 
to Noether's theorem and the respective manifestations of the symmetries in nature. 

The SUy(3) and Uy(l) symmetries ensure the conservation of isospin and baryon number, 
respectively, while the SUa(3) and Ua(l) symmetries are transformations that involve the 
75 matrix and therefore alter the parity of the state in which they operate. For the sake of 
uniformity throughout the text we will designate by chiral symmetry the SU^(3) symmetry 
and by axial symmetry the 11^(1) symmetry. 

From the experimental point of view the manifestation of chiral symmetry would be the 
existence of parity doublets, i.e., a multiplet of particles with the same mass and opposite 
parity for each multiplet of isospin (the chiral partners), in the hadronic spectrum; that 
situation is not verified. Similarly, if the symmetry U^(l) would manifest itself, the existence 
of a partner with opposite parity to each hadron should be observed experimentally. As 
neither of these situations is observed in the hadron spectrum, these symmetries must be 
somehow broken. 

Concerning the SU^4(3) symmetry, the theory must contain a mechanism for the sponta- 
neous breaking of chiral symmetry, which represents a transition to an asymmetric phase. 
This is closely related to the existence of non-zero quark condensates, (qq), which are not 
invariant under SU^(3) transformations and therefore act as order parameters for the spon- 
taneously broken chiral symmetry. 

According to the Goldstone theorem, the spontaneous breaking of a continuous global 
symmetry implies the existence of a particle with zero mass, the Goldstone boson. In the 



case we are considering, the symmetry breaking is closely related to the appearance of 
eight degenerate Goldstone bosons with zero mass. As matter of fact the pions where the 
first mesons associated to Goldstone bosons due to their small mass. Indeed, if compared 
to the mass of the nucleon one has M n /M^ = 0.15. To reproduce the meson spectrum 
it is also necessary that the theory incorporates a mechanism which explicitly breaks the 
chiral symmetry: the Lagrangian must include perturbative terms that break ab initio the 
symmetry allowing the lifting of the degeneracy in the pseudoscalar mesons spectrum. These 
terms are the current quark masses m u m^ and m s . 

Nevertheless, the fundamentals of this process are not yet completely understood at the 
level of QCD. However, the chiral symmetry breaking, the generation of the constituent 
quark masses and the concomitant appearance of Goldstone bosons are explained in a very 
satisfactory way in certain theories and physical models like NJL type models. 

Another symmetry effect occurs in the 77 — 77' sector. The pseudoscalar 77' meson is 
not of the Goldstone type. Classically, the QCD Lagrangian with Nf = 3 massless flavors 
possesses the Ul(3) xU^(3) symmetry. In fact, the Adler-Jackiw-Bell Ua(1) anomaly breaks 
this symmetry to SUl(3)x SUr(3) and yields a large mass for the 77' even if quark masses 
are zero. The origin of this anomaly are the instantons in QCD as explained in the next 
subsection. 

In conclusion, QCD with three massless quarks possesses the SUl(3)x SU#(3) global 
chiral symmetry. At low energy, the pertinent degrees of freedom are mesons indicating that 
QCD experiences a confined/deconfined phase transition when energy increases. The study 
of mesons give us insight on the non-perturbative vacuum of QCD at low energy. 

C. 11^(1) Symmetry Breaking 

Now we will analyze the Ua(1) symmetry. It is known that in the chiral limit, and at the 
classical level, Cqcd is also invariant under the axial transformation (see Tabled]). 

As already mentioned, the absence in nature of chiral partners with the same mass, related 
with the U^(l) symmetry, opens the possibility that U^(l) symmetry is also spontaneously 
broken, similarly to what happens with the SILa(3) symmetry. Consequently, there should 
exist another pseudoscalar Goldstone boson. S. Weinberg estimated the mass of this particle, 



outside the chiral limit, at about y/^M^ 60]. 
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From all hadrons known in nature, the only candidates with the correct quantum num- 
bers were the 7/(549) and 7/(985). However, both violate the Weinberg estimate because 
their masses are too high: the r/(985) has a mass comparable to the mass of nucleons ("77' 
problem") and, furthermore, the 7/(549) is identified as belonging to the octet of pseudoscalar 
mesons. If a particle with the characteristics pointed out by Weinberg does not exist, then 
where is the ninth Goldstone boson? If this Goldstone boson does not exist there is no 
spontaneous breaking of Ua(1) symmetry. 

In 1976 G. 't Hooft suggested that the Ua(1) symmetry does not exist at the quantum 



level, being explicitly broken by the axial anomaly 60[ which can be described at the semi 



classical level by instantons [6lLl62|. Instantons are AD topologically stable soliton solutions 
of QCD. Other well known solitons are 2D Abrikosov vortices in superconducting material or 
3-D magnetic monopole in QCD. The instantons "transform" left handed fermions into right 
handed ones (and conversely). The transition induced by this solution of QCD has a non 
zero axial charge variation: AQ 5 = ±2Nf, i.e., one left handed fermion (+1 in axial charge) 
transforms in a right handed fermion (—1). To mimic this interaction in a purely fermionic 
effective theory, 't Hooft proposed the following interaction Lagrangian, the so-called '"t 
Hooft determinant": 

dust = Ke l0mat det (i) R (x)i) L (x)) + h.c. (5) 

flavor 

This interaction "absorbs" Nf left helicity fermions and converts them to right handed ones 
(and conversely). In the following, we will take 8i nst = 0. 

In this context, as suggested by 't Hooft, the instantons can play a crucial role in breaking 
explicitly the Ua(1) symmetry giving to the 77' a mass of about 1 GeV outside the chiral 
limit. This implies that the mass of 77' has a different origin than the other masses of the 
pseudoscalar mesons and can not be seen as the missing Goldstone boson due to spontaneous 
breaking of Ua(1) symmetry. Consequently, the Ua(1) anomaly is very important once it is 
responsible for the flavor mixing effect that removes the degeneracy among several mesons. 

Finally, an interesting open question also related to the symmetry that we wish to address 
is whether both chiral S\J(Nf)xS\J(Nf) and axial Ua(1) symmetries are restored in the 
high temperature/density phase and which observables could carry information about these 
restorations. 



D. The Polyakov Loop and the Z3 Symmetry Breaking: Pure Gauge Sector 



In this Section, following the arguments given in [38|, |39[ , we discuss how the deconfine- 
ment phase transition in a pure SU(iV c ) gauge theory can be conveniently described through 
the introduction of an effective potential for the complex Polyakov loop field, which we define 
in the following. 

Since we want to study the SU(iV c ) phase structure, first of all an appropriate order 
parameter has to be defined. For this purpose the Polyakov line 



L (x) = V exp 



f 
1 dr v4 4 (2, r) 

Jo 



(6) 



is introduced. In the above, A4 = iA is the temporal component of the Euclidean gauge 
field (A, Ai), in which the strong coupling constant g has been absorbed, V denotes path 
ordering and the usual notation (3 = 1/T has been introduced with the Boltzmann constant 
set to one (ks = 1). 

The Polyakov line L (x) can be described as an operator of parallel transport of the gauge 
field A 4 (x, r) into the direction r. One way to understand why this quantity is indeed a 
parameter that can distinguish between confined or deconfined phase is to consider the two 
extremal behavior of L: a color field A4 at the point (x, r) will be transformed into the color 
field L x A4 after being transported into the direction r. If L — >• 1 it means that nothing 
affects the propagation of the field: the medium is in its deconfined phase. At the contrary 
L — > indicates that the color field cannot propagate in the medium: it is confined. 

Another way to see this point is to consider the variation of free energy when an infinitely 
massive (hence static) quark (that acts as a test color charge) is added to the system. To 
this purpose let us introduce the Polyakov loop. When the theory is regularized on the 
lattice, it reads: 

/(£) = ^-TrL(f) (7) 

and it is a color singlet under SU(iV c ), but transforms non-trivially, like a field of charge 
one, under Z^. Its thermal expectation value can then be chosen as an order parameter 



for the deconfinement phase transition 63H65J : in the usual physical interpretation 66l. l67| 



(l(x)) is related to the change of free energy occurring when a heavy color source in the 
fundamental representation is added to the system. One has: 

(/(f)) = e -0AF o (a) (8) 
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In the Zjv c symmetric phase, (l{x)) = 0, implying that an infinite amount of free energy is 
required to add an isolated heavy quark to the system: in this phase color is confined. 
In the case of the SU(3) gauge theory, the Polyakov line L{x) gets replaced by its gauge 

Note that ((L{x))) 



covariant average over a finite region of space, denoted as ((L(x))) 38 
in general is not a SU(iV c ) matrix. The Polyakov loop field 

1 



$(x) 



((Km 



N, 



Tr c «L(x)» 



(9) 



is then introduced. 

Following the Landau-Ginzburg approach, a Z 3 symmetric effective potential is defined 
for the (complex) $ field, which is conveniently chosen to reproduce, at the mean field level, 



results obtained in lattice calculations 



38 



39 



43| . In this approximation one simply sets 



the Polyakov loop field $(x) equal to its expectation value $ = const., which minimizes the 
potential. 



Concerning the effective potential for the 
in the literatur e | 68l -l70l] ; the one proposed in 



results 



69 



71 



72 



^complex) field $, different choices are available 
see Equation ( TT0]) ) is known to give sensible 



and will be adopted in our parametrization of the PNJL model that will 
be presented in Section IIIII In particular, this potential reproduces, at the mean field level, 
results obtained in lattice calculations as it will be shown. The potential reads: 

W($,$;T) a(T) 



rp4 



-$$ + 6(T)ln[l - 6$$ + 4(<T + $ 3 ) - 3($$)' 



where 



a(T) = a + a 1 I -£j + a 2 l-^J and b(T) = b 3 (-£ 



(10) 



(11) 



The effective potential exhibits the feature of a phase transition from color confinement 
(T < To, the minimum of the effective potential being at $ = 0) to color deconfinement 
(T > Tq, the minima of the effective potential occurring at $ ^ 0). 

The parameters of the effective potential U are given in Table [Til These parameters have 
been fixed in order to reproduce the lattice data for the expectation value of the Polyakov 
loop and QCD thermodynamics in the pure gauge sector [73|, |74| . 

The parameter To is the critical temperature for the deconfinement phase transition within 
a pure gauge approach: it was fixed to 270 MeV, according to lattice findings. Different 
criteria for fixing T may be found in the literature, like in 75|, where an explicit Nf 
dependence of T is presented by using renormalization group arguments. However, it is 
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a 


ai 


a2 


&3 


3.51 


-2.47 


15.2 


-1.75 



TABLE II. Parameters for the effective potential in the pure gauge sector (Equation (|10p ). 



noticed that the Polyakov loop computed on the lattice with(2-r-l) flavors and with fairly 



realistic quark masses is very similar to the SU/(2) case |74j. Hence, it was chosen to 
keep for the effective potential U ($, $; T\ the same parameters which were used in SU/(2) 
PNJL 69[, including To = 270 MeV. The latter choice ensures an almost exact coincidence 
between chiral crossover and deconfinement at zero chemical potential, as observed in lattice 
calculations. We notice, however, as it will be seen in Section IIVt that a rescaling of T to 
210 MeV may be needed in some cases in order to get agreement between model calculations 
and thermodynamic quantities obtained on the lattice. Let us stress that this modification 
of To (the only free parameter of the Polyakov loop once the effective potential is fixed) is 
essentially done for rescaling reasons (an absolute temperature scale has not a very strong 
meaning in these kind of Ginzburg-Landau model based on symmetry) but does not change 
drastically the physics. 

With the present choice of parameters, $ and $ are always lower than 1 in the pure 
gauge sector. In any case, in the range of applicability of our model, T < (2 — 3)T C , there is 
a good agreement between the results and the lattice data for $. 

To conclude this section, some comments are in order on the number of parameters and 
the nature of the Polyakov field. 

It should be noticed that the NJL parameters and the Polyakov potential ones are not 
on the same footing. Whereas the NJL parameters are directly related in a one to one 
correspondence with a physical quantity, the Polyakov loop potential is there to insure 
that the pure gauge lattice expectations are recovered. Hence the potential for the loop 
can be viewed as a unique but functional parameter. The details of this function are not 
very important to study the thermodynamic as soon as the potential reproduces the lattice 
results. In order to clarify this point, we remember that in a previous calculations in the 
Nf = 2 case 76[ we used two kind of potentials and obtained very small differences in what 
concerns thermodynamics; however, if one calculates susceptibilities with respect to $ the 
log potential used here has to be preferred in order not to have unphysical results. The only 
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true parameter is the pure gauge critical temperature To that fixes the temperature scale of 
the system. However we would like to point out that it is expected in the Landau-Ginzburg 
framework that a characteristic temperature for a phase transition will not be a prediction: 
one needs to fix the correct energy scale somehow and it is the role of this parameter. Hence 
we will allow ourselves to change this parameter in the calculation of some observables in 
order to compare our results with lattice QCD expectations. 

Finally we want to stress that, contrarily to the full Landau-Ginzburg effective field 
approach, the Polyakov loop effective field is not a dynamical degree of freedom due to the 
the lack of dynamical term in the Polyakov loop potential. Hence it is a background gauge 
field in which quarks will propagate. Anyway the potential mimics a pressure term that has 
the correct magnitude and temperature behavior in order to get the Stefan-Boltzman limit 
for the gluonic degrees of freedom and that explains the success of the model. 

E. QCD Phase Diagram 

Understanding the QCD phase structure is one of the most important topics in the physics 
of strong interactions. The developed effort on both the theoretical point of view — using 
effective models and lattice calculations — and the experimental point of view is one of the 
main goals of the heavy ion collisions program; it has proved very fruitful, shedding light on 
properties of matter at finite temperatures and densities. 

From the theoretical point of view, the physics of the transition from hadronic matter to 
the QGP is well established for baryonic chemical potential \xb = 0. Latest lattice results for 
(2+1) flavors show a crossover transition at T c = (170 — 200) MeV (using m q = 0.05 m s ) J77J. 
Near the critical temperature, T c , the energy density (and other thermodynamic quantities) 
show a strong growth, signaling the transition from a hadronic resonance gas to a matter 
of deconfined quarks and gluons. As matter of fact, the rapid rise of the energy density 
is usually interpreted to be due to deconfinement i.e., liberation of many new degrees of 
freedom. In the limit where quark masses are infinitely heavy the order parameter for the 
deconfinement phase transition is the Polyakov loop. A rapid change in this quantity is also 
an indication for deconfinement even in the presence of light quarks (77| . 

Concerning chiral symmetry, and considering the chiral limit, it is expected a chiral 
transition with the corresponding order parameter being the quark condensate: the quark 
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condensate vanishes at the critical temperature T c and a genuine phase transition takes 
place. Even away from the chiral limit, when the quark masses are finite, it is expected, in 
the transition region, a "crossover" where the quark condensates rapidly drop indicating a 
partial restoration of the chiral symmetry J77l |. 

This rises the interesting question whether the restoration of chiral symmetry and the 
transition to the QGP occurs at same time. In [78| it was proposed that the evidence of 
restoration of chiral symmetry is a sufficient condition to demonstrate the existence of a new 
state of matter, but it is not a necessary condition for the discovery of the QGP. In fact, 
as already pointed out, most of lattice results shows a tendency for the restoration of chiral 
symmetry to happen simultaneously with deconfinement, but this issue is not definitively 
determined from a theoretical point of view. 

At finite temperature and chemical potential the most common three-flavor phase diagram 
shows a first order boundary of the chiral phase transition separating the hadronic and quark 
phases. This first-order line starts at non-zero chemical potential and zero temperature and 
will finish in a point, the critical end point (T CEP , fi^ EP ), where the phase transition is of 
second order. As the temperature increases and the chemical potential decreases, the QCD 
phase transition becomes a crossover. 

Although recent lattice QCD results by de Forcrand and Philipsen question the existence 
of the CEP 79M81J. this critical point of QCD, proposed at the end of the eighties 45l-l48| . 



is still a very important subject of discussion nowadays 51|. The search of the QCD CEP is 



one of the main goals in "Super Proton Synchrotron" (SPS) at CERN [82| and in the next 
phase of the "Relativistic Heavy Ion Collider" (RHIC) running at BNL. 

When matter at high density and low temperature is considered, it is expected that this 
type of matter is a color superconductor where pairs of quarks condense ( "diquark conden- 
sate"): this color superconductor is a degenerate Fermi gas of quarks with a condensate 
of Cooper pairs near the Fermi surface that induces color Meissner effects. With the color 
superconductivity it was realized that the consideration of diquark condensates opened the 
possibility for a wealth of new phases in the QCD phase diagram. At the highest densi- 
ties/chemical potentials, where the QCD coupling is weak, the ground state is a particularly 
symmetric state, the color-flavor locked (CFL) phase: up, down, and strange quarks are 
paired in a so-called color-flavor locked condensate. At lower densities the CFL phase may 
be disfavored in comparison with alternative new phases which may break translation and/or 
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rotation invariance 36l. |83|. 



Recently it was argued that some features of hadron production in relativistic nuclear 
collisions, mainly at SPS-CERN energies, can be explained by the existence of a new form of 
matter beyond the hadronic matter and the QGP: the so-called quarkyonic matter 
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861 ] . It was also suggested that these different types of matter meet at a triple point in the 
QCD phase diagram: both the hadronic matter, the QGP, and the quarkyonic matter all 
coexist [87 1 . 



From the experimental point of view, different parts of the QCD phase diagram have been 
investigated for different beam energies. For high energies, corresponding to the energies of 
RHIC-BNL and LHC-CERN ("Large Hadron Collider"), the matter is produced at high 



temperature and low baryonic density 78 



On the other hand, intermediate energy, 



corresponding to the energies of "Alternating Gradient Synchrotron" (AGS) in BNL or 
in SPS-CERN 82|, can be studied through matter formed at moderate temperature and 
baryonic density. Finally, the energies achieved in the GSI Heavy Ion Synchrotron (SIS) 
90| and "Nuclotron-based Ion Collider fAcility" (NICA) at JINR [9l| will allow to obtain 
matter at low temperatures and high baryonic densities. For example the main goal of 
the research program on nucleus-nucleus collisions at GSI is the investigation of highly 
compressed nuclear matter. Matter in this form exists in neutron stars and in the core of 
supernova explosions. 

In Figure [1] are indicated the points where the hadrons cease to interact with each others 
(freeze-out points) for the different beam energies. These points are related to final state of 
the expansion of the fireball. 

It should also be noted that atomic nuclei by themselves represent a system at finite 
density and zero temperature. At normal nuclear density it is estimated that the quark 
condensate undergoes a reduction of 30% [93| so that the resulting effects in the chiral order 
parameter can be measured in beams of hadrons, electrons and photons in nuclear targets. 

If the evidence for deconfinement can be found experimentally, then the search for man- 
ifestations of chiral symmetry restoration will be, in the near future, one of the main goals 
once the investigation of the properties of matter can provide a clear evidence for changes 



in the fundamental vacuum of QCD, with far-reaching consequences 78 
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FIG. 1. The phase diagram from the experimental point of view: freeze-out points for the different 
beam energies are showed 921 ] . 



III. THE PNJL MODEL WITH THREE FLAVORS 

The Polyakov-Nambu-Jona-Lasinio (PNJL model) that we want to build is an effective 
model for QCD, written in term of quark degrees of freedom. In the next paragraphs we 
will present the ingredients we need to build our model namely: a massive Dirac Lagrangian 
together with a four quark chirally invariant interaction (the original NJL Lagrangian with 
a small mass term that breaks explicitly the chiral symmetry, has observed in the mass 
spectrum); the so-called 't Hooft interaction that reproduces the interaction of quarks with 
instantons and finally a Polyakov loop potential that mimic the effect of the pure gauge 
(Yang-Mills) sector of QCD on the quarks. 

A. Nambu Jona-Lasinio Model with Anomaly and Explicit Symmetry Breaking 

Phase transitions are usually characterized by large correlation lengths, i.e., much larger 
than the average distance between the elementary degrees of freedom of the system. Effective 
field theories then turn out to be a useful tool to describe a system near a phase transition. 
In particular, in the usual Landau- Ginzburg approach, the order parameter is viewed as 
a field variable and for the latter an effective potential is built, respecting the symmetries 
of the original Lagrangian. The existence of a phase transition between two sectors where 
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the chiral symmetry is spontaneously broken or restored (transition associated to the quark 



condensate that acts as an order parameter) and t 



re Ginzburg-Landau theory suggests to 



use the symmetry motivated NJL Lagrangian [35|, [36|, [59|, |94( for the description of the 
coupling between quarks and the chiral condensate in the scalar-pseudoscalar sector. The 
associated Lagrangian which complies with the underlying symmetries of QCD described in 
the previous section reads: 



£>NJL — Cq + Ci + Cq (12) 

A) = q(iYd^ - m)q (13) 

1 8 

£ 4=2^ £[UA a ?) 2 + (^7 5 A a g) 2 ] (14) 

a=0 

C 6 = g D {det [q (1 + 75) q] + det [q (1 - 75) <?]} (15) 

In the above q = (u, d, s) is the quark field with three flavors (TV/ = 3) and three colors 
(N c — 3), rh — dia,g(m u ,md,m s ) is the current quark mass matrix, and A a are the flavor 
811/(3) Gell-Mann matrices (a = 0, 1, . . . , 8), with A = J\ I. 

As it is well known, this Lagrangian is invariant under a global — and nonlocal — color 
symmetry SU(N C = 3) and lacks the confinement feature of QCD. 

£4 has Ul(3) x Ur(3) chiral symmetry and Cq is a six quarks interaction, the 't Hooft 
determinant to reproduce {7a (1) anomaly. The £7^(3) x Ur{3) symmetry of £ + £ 4 is broken 
to SU L {3) x SU R {3) by £ 6 . 

The current quark mass matrix rh is in general non degenerate and explicitly breaks 
chiral symmetry SUl(3) x SUr(3) to S77/(3) or its subgroup. In the following we will take 
m u = rrid hence keeping exact the isospin symmetry. 

Let us notice that NJL has no build-in confinement and is non-renormalizable thus it 
requires the introduction of a cutoff parameter A. 



The parameters, obtained by following the methodology of 95| . are: m u = m^ = 5.5 
MeV, m s = 140.7 MeV, g s A 2 = 3.67, g D A 5 = -12.36 and A = 602.3 MeV, which are fixed 
to reproduce the values of the coupling constant of the pion, /„. = 92.4 MeV, and the meson 
masses of the pion, the kaon, the 77 and 77', respectively, M n = 135 MeV, M K = 497.7 MeV, 



M v = 514.8 MeV and M v , = 960.8 MeV. 
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B. Coupling between Quarks and the Gauge Sector: The PNJL Model 

In the presence of dynamical quarks the Z3 symmetry is explicitly broken: the expectation 
value of the Polyakov loop still serves as an indicator for the crossover between the phase 
where color confinement occurs ($ — > 0) and the one where color is deconfined ($ — > 1). 

The PNJL model attempts to describe in a simple way the two characteristic phenomena 
of QCD, namely deconfinement and chiral symmetry breaking. 

In order to describe the coupling of quarks to the chiral condensate, we start from a 
NJL description of quarks; they are coupled in a minimal way to the Polyakov loop, via a 
covariant derivative. Hence our calculations are performed in the framework of an extended 
SU/(3) PNJL Lagrangian, which includes the 't Hooft instanton induced interaction term 



that breaks the Ua(1) symmetry and 



temporal background gauge field $ 
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;he quarks are coupled to the (spatially constant) 
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97j . The Lagrangian reads: 



1 8 

Cpnjl =q{i^D^-m)q+-g s ^[(<?A a g) 2 + {qil,\ a qf} 

a=0 

+ ^{det[g(l + 7 5)g] + det[g(l- 75 )g]}-W($[A],<l[A];T) (16) 

The covariant derivative is defined as D yi = d^ — iA 11 , with A^ = SqAq (Polyakov gauge); in 
Euclidean notation Aq = —iA^. The strong coupling constant g is absorbed in the definition 
of A^(x) = gA^(x)^f, where A% is the (SU C (3)) gauge field and A a are the (color) Gell-Mann 
matrices. 

At T = 0, it can be shown that the minimization of the grand potential leads to $ = $ = 
0. So, the quark sector decouples from the gauge one, and the model is fixed as referred in 
the previous subsection. 

Some remarks are in order concerning the applicability of the PNJL model. It should 
be noticed that in this model, beyond the chiral point-like coupling between quarks, the 
gluon dynamics is reduced to a simple static background field representing the Polyakov 



loop (see details in [44, l68|). This scenario is expected to work only within a limited range 
of temperatures, since at large temperatures transverse gluons are expected to start to be 
thermodynamically active degrees of freedom and they are not taken into account in the 
PNJL model. We can assume that the range of applicability of the model is roughly limited 



to T < (2— 3)T C , since, as concluded in j98|, transverse gluons start to contribute significantly 



for T > 2.5 T c , where T c is the deconfinement temperature. 
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C. Grand Potential in the Mean Field Approximation 

Our model of strongly interacting matter can simulate either a region in the interior of 
neutron stars or a dense fireball created in a heavy-ion collision. In the present work, we 
focus our attention on the last type of systems. Bearing in mind that in a relativistic heavy- 
ion collision of duration of about 10 -22 s, thermal equilibration is possible only for processes 
mediated by the strong interaction rather than the full electroweak equilibrium, we impose 
the condition fi e = and use the chemical equilibrium condition \x u = fid — Us — /•*• This 
choice allows for isospin symmetry, m u = m d and approximates the physical conditions at 
RHIC. 

We will work in the (Hartree) mean field approximation. In this context, the quarks can 
be seen as free particles whose bare current masses rrii are replaced by the constituent (or 
dressed) masses M,. 

The quark propagator in the constant background field A4 is then: 

Si(p) =-(rf-M i + 7o (/i - iM))- 1 (17) 

:or a fermion. 



35l | ) to obtain ef- 



In the above, po = i^n an d uj n = (2n + l)irT is the Matsubara frequency 

Within the mean field approximation, it is straightforward (see Ref. |£ 
fective quark masses from the Lagrangian (1161) ; these masses are given by the so-called gap 
equations: 

M i = m l - 2g s (£$) - 2g D {q^) (q k q k ) (18) 

where the quark condensates (ftft), with i,j, k = u,d,s (to be fixed in cyclic order), have to 
be determined in a self-consistent way as explained in the Appendix. 

The PNJL grand canonical potential density in the SU/(3) sector can be written as 

n = Q^,^,M i ;T, f i)=U(^,^,T)+g s ^ (q^f + Ag D {q u q u ) (q d q d ) (q s q s ) 

{i=u,d,s} 

- 2N < E / a t^- 2T E /^(«+«) (19) 



where E { is the quasi-particle energy for the quark i: E{ = ^p 2 + M 2 , and z$ and z^ are 
the partition function densities. 
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The explicit expression of z$ and z$ are given by: 

4(£,) = Tr c ln [l + l)e~P E A = In {l + 3 ($ + Se^** ) e"^ + e -3 "** } (20) 
z*(£i) = Tr c ln 1 + Le~^\ = In {l + 3 ($ + $e^ + ) e"^ + + e" 3 ^} (21) 

where J5J = ^ =)= /i, the upper sign applying for fermions and the lower sign for anti- 
fermions. The technical details are postponed to the Appendix but as a conclusion on this 
topic, a word is in order to describe the role of the Polyakov loop in the present model. 
Almost all physical consequences of the coupling of quarks to the background gauge field 
stem from the fact that in the expression of z$, $ or $ appear only as a factor of the one- or 
two-quarks (or antiquarks) Boltzmann factor, for example e~^ Ep and e~ 2l3Ep . Hence when 
$, $ — > (signaling what we designate as the "confined phase") only e~ 3 @ Ep remains in the 
expression of the grand potential, leading to a thermal bath with a small quark density. At 
the contrary $, $ — > 1 (in the "deconfined phase") gives a thermal bath with all 1-, 2- and 
3-particle contributions and a significant quark density. Other consequences of this coupling 
of the loop to the 1- and 2-particles Boltzmann factor will be discussed in Section IVIIII 

D. Equations of State and Response Functions 

The equations of state can be derived from the thermodynamical potential £l(T, ft). This 
allows for the comparison of some of the results with observables that have become accessible 
in lattice QCD at non-zero chemical potential. 



As usual, the pressure p is defined such as its value is zero in the vacuum state 36[ and, 
since the system is uniform, we have 

p(T,Ai) = -^[n(T,Ai)-n(0,0)] (22) 

where V is the volume of the system. 

The relevant observables are the baryonic density 

'^"HfrX (23) 

and the (scaled) "pressure difference" given by 

Ap(T,/x) = p(T,/x)-p(T,0) 

rp A rp A \ ) 
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Due to the relevance for the study of the thermodynamics of matter created in relativistic 
heavy-ion collisions, it is interesting to perform an analysis of the isentropic trajectories. The 
equation of state for the entropy density, s, is given by 

«CT,A0 = (H) (25) 

and the energy density, e, comes from the following fundamental relation of thermodynamics 

e(T,p) = Ts{T,p) + pp B {T,p) - p(T, p) (26) 

The energy density and the pressure are defined such that their values are zero in the 



vacuum state |36| 



The baryon number susceptibility, xb, and the specific heat, C, are the response of 
the baryon number density, pb{T,p), and the entropy density, s(T,p), to an infinitesimal 
variation of the quark chemical potential p and temperature, given respectively by: 

*■-(*), - '-£(©. 

These second order derivatives of the pressure are relevant quantities to discuss phase 
transitions, mainly the second order ones. 



E. Model Parameters and Regular izat ion Procedure 

As already seen, the pure NJL sector involves five parameters: the coupling constants 
gs and gp, the current quark masses m u = m^, rn s and the cutoff A defined in Section UTT1 
These parameters are determined in the vacuum by fitting the experimental values of several 
physical quantities. We notice that the coupling constant gs and go and the parameter A are 
correlated. For instance, if we increase gs in order to provide a more significant attraction 
between quarks, we must also increase the cutoff A in order to insure a good agreement with 
experimental results. In addition, the value of the cutoff itself does have some impact as far 
as the medium effects in the limit T = are concerned. Also the strength of the coupling 
go has a relevant role in the location of the CEP, as it will be discussed. 

In fact, the choice of the parametrization may give rise to different physical scenarios at 
T = and p B = p ^ [36j , even if they give reasonable fits to hadronic vacuum observables 
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and predict a first order phase transition. The set of parameters we used insures the stability 
conditions and, consequently, the compatibility with thermodynamic expectations. 

On the other hand, the regularization procedure, as soon as the temperature effects are 
considered, has relevant consequences on the behavior of physical observables, namely on 
the chiral condensates and the meson masses 99|. Advantages and drawbacks of t hese 



regularization procedures have been discussed within the NJL 99[ and PNJL 76|, Il00 | 
models. We remind that one of the drawbacks of the regularization that consists in putting 
a cutoff only on the divergent integrals is that, at high temperature, there is a too fast 
decrease of the quark masses that become lower than their current values. This leads to 
a non physical behavior of the quark condensates that, after vanishing at the temperature 
(T e ff) where constituent and current quark masses coincide, change sign and acquire non- 
zero values again. To avoid this unphysical effects we use the approximation of imposing by 
hand the condition that, above T e ff, Mi = rrii and {qtqi) = 0. 

IV. EQUATION OF STATE AT FINITE TEMPERATURE 
A. Characteristic Temperatures 

At zero temperature and chemical potential, the chiral symmetry of QCD is explicitly 
broken. It is expected that chiral symmetry will be restored at high temperature; hence a 
phase transition occurs separating the low and the high temperature regions. This phenom- 
ena may be realized in high energy heavy ion collision experiments. 

In this case (fiB — and T^O), where a crossover transition occurs in the PNJL model, a 
unique critical temperature cannot be defined, but one of many characteristic temperatures 
for the transition may be used. Of course, these temperatures should coincide in the limit 
where the transition becomes of second order (in the chiral limit, for example, when one is 
concerned with the restoration of chiral symmetry) . 

We start our analysis by identifying the characteristic temperatures which separate the 



different thermodynamic phases in the PNJL model [44j , using the regularization that allows 
high momentum quark states at finite temperature (the cutoff is used only in the integrals 
that are divergent and A — > oo in the ones that are already convergent because of the 
thermal distributions). 
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T [MeV] 


T x [MeV] 


If [MeV] 


T c [MeV] 


270 


222 


210 


216 


210 


203 


171 


187 



TABLE III. Characteristic temperatures with two different values for Tg. 



Let us analyze the behavior of the quark masses and the field $. The characteristic tem- 
perature related to the deconfinement phase transition is T* and the chiral phase transition 
characteristic temperature, T x , signals partial restoration of chiral symmetry. These tem- 
peratures are chosen to be, respectively, the inflexion points of the "quasi" order parameter 
$ and of the chiral condensate (q u q u )] as in 68J, we d efine T c as the average of the two 
transition temperatures, Tf and T x . As shown in 100J], the present regularization lowers 
the characteristic temperatures and decreases the gap T* — T x , leading therefore to better 
agreement with lattice results. 

Let us remark that this section is also devoted to the study of thermodynamic quantities 
(pressure, energy per particle and entropy) for which a rescaling of the temperature T is 
needed, in order to get a better agreement with lattice results for these quantities. Therefore, 
following the argumentation presented in 681 ]. here we will use the reduced temperature T c 
by rescaling the parameter To from 270 to 210 MeV (let us stress that this rescaling is only 
done for the remainder of this section). Results for the characteristic temperatures with 
T = 210 MeV and T = 270 MeV are shown in Table ED 

It should be noticed that, when T = 210 MeV is used, we loose the almost perfect 
coincidence of the chiral and deconfinement transitions (they are shifted relative to each 
other by about 32 MeV) and we have T c = 187 MeV within the range expected from lattice 
calculations 771 ] . However, the behavior of the relevant physical quantities is qualitatively 
the same whether T = 270 MeV or T = 210 MeV. 

An interesting point to be noticed in Figure [21 lower panel, is that the inflexion points 
of the strange and non-strange quark condensates are close and slightly separated from the 
inflexion point of $, although the condensates have small peaks around the inflexion point 
of $. The latter effect is r elate d to the fastening of the phase transition induced by the 
Polyakov loop as shown in 1021 ] (it obviously does not signal a phase transition since the 
variation of the order parameter around this temperature is small). The Polyakov loop leads 
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FIG. 2. Upper part: quark masses (left panel) and quark condensates (right panel) in PNJL 
model as functions of the temperature; the Polyakov loop field $ is also shown. At T e ff = 345 
MeV, M{ = 77ij. Lower part: derivatives of the quark condensates and of the Polyakov loop field 
<3? (left panel); the topological susceptib ility , x> (right panel) in the PNJL model, compared to 



corresponding lattice results taken from 



1011 ] . 



to a faster decrease of the quark masses around T c and the present regularization enhances 
this effect, even at temperatures higher than T c . At T e ff = 345 MeV the regularization 
is responsible for the full restoration of the chiral symmetry that was dynamically broken: 
the quark masses go to their current values and the quark condensates vanish (see Figure 
El left and right panels) . As already shown in the framework of the pure NJL model 99] 



and in the PNJL model 



100J, the effect of allowing high momentum quark states is stronger 



for the strange quark mass. Indeed, with the conventional regularization, the non-strange 
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constituent quark mass at high temperature is already very close to its current mass, the new 
regularization only enhancing this behavior. Differently, the strange quark mass is always 
far from its current value, unless we allow high momentum quarks to be present and, in that 
case, its constitu ent mass decreases very substantially and comes to its current value. 



As shown in 100], at T e ff the behavior of some given observables signals the effective 
restoration of chiral and axial symmetry: the masses of the meson partners of both chiral 
and axial symmetry are degenerated and the topological susceptibility vanishes as we can 
see from Figure |2] (lower right panel). 

B. Thermodynamic Quantities 

In the limit of vanishing quark chemical potential, significant information on the phase 
structure of QCD at high temperature is obtained from lattice calculations. The transition 
to the phase characteristic of this regime is related with chiral and deconfinement transitions 
which are the main features of our model calculation. 

In Figure [31 we plot the scaled pressure, the energy and the entropy as functions of the 

temperature compared with recent lattice results (see Reference [77|). Since the transition 

to the high temperature phase is a rapid crossover rather than a phase transition, the 

pressure, the entropy and the energy densities are continuous functions of the temperature. 

We observe a similar behavior in the three curves: a sharp increase in the vicinity of the 

transition temperature and then a tendency to saturate at the corresponding ideal gas limit. 

Asymptotically, the QCD pressure for Nf massless quarks and (N% — 1) massless gluons is 

given (/i B = 0) by: 

?*» = (NZ-l)- + N e N,—, (28) 

T 4 v c ; 45 ' 180 ' K J 

where the first term denotes the gluonic contribution and the second term the fermionic one. 

The results follow the expected tendency and go to the free gas values (or Stefan- 

Boltzmann limit), a feature that was also foun d with this type of regularization in the 



context of the SU(2) PNJL model 68|, Il03l . Il04j . The inclusion of the Polyakov loop effec- 



tive potential W($, $; T) (it can be seen as an effective pressure term mimicking the gluonic 
degrees of freedom of QCD) is required to get the correct limit (indeed in the NJL model 
the ideal gas limit is far to be reached due to the lack of gluonic degrees of freedom) . 

The inclusion of the Polyakov loop and the regularization procedure are essential to obtain 
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FIG. 3. Scaled pressure (p), energy per particle (e), and entro py (s) as a function of the temperature 
at zero chemical potential. The data points are taken from [77( ■ The pressure reaches 66% of the 
strength of the Stefan-Boltzmann value at T = 300 MeV, a value which attains 85% at T = 400 
MeV. 



the required increase of extensive ther mody namic quantities, insuring the convergence to the 
Stefan-Boltzmann (SB) limit of QCD |105| . Some comments are in order concerning the role 
of the regularization procedure for T > T c . In this temperature range, due to the presence of 
high momentum quark states, the physical situation is dominated by the significant decrease 
of the constituent quark masses by the qq interactions. This allows for an ideal gas behavior 
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of almost massless quarks with the correct number of degrees of freedom. 

Let us notice that, just below T c , the pressure and the energy fail to reproduce the lattice 
points: for example there is a small underestimation of the pressure and energy in the 
model calculations. It is known that the lack of mesonic correlations in the PNJL model 



is responsible for, at least, a fraction of this discrepancy. As matter of fact, in [106J ] the 
authors have evaluated the pion and sigma contributions to the pressure by calculating the 
ring sum in the SU(2) nonlocal PNJL model. They have shown that the pionic contribution 
dominates the pressure in the low temperature region, shifting the respective curve to higher 
values. It is expected that the introduction of mesonic correlations in the SU(3) PNJL model 
has the same effect for the pressure. 

V. PHASE DIAGRAM AND THE LOCATION OF THE CRITICAL END POINT 

Although recent lattice QCD results by de Forcrand and Philipsen question the existence 



of the CEP 
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8l|, this critical point of QCD, proposed at the end of the eighties [4a-l48j . 
is still a very important subject of discussion nowadays 51] . 

We remember that the TCP separates the second order transition at high temperature 
and low chemical potential, from the first order transition at high chemical potential and 
low temperature. If the second order transition is replaced by a smooth crossover, a CEP 
which separates the two lines is found. In order to determine and elucidate the nature of 
the phase transition the relevant thermodynamic quantities are studied, starting with zero 
temperature and finite chemical potential. 

A. Phase Transition at Zero Temperature 

Figure 0] illustrates the properties of cold quark matter as predicted by the model cal- 
culation. In the left panel the negative of the thermodynamic potential is plotted as a 
function of the chemical potential, showing the presence of branches with stable, metastable 
and unstable solutions. This behavior is obtained in the domain \j? B < jjlb < [a b , where 
the gap equations ( fl~8l) have three solutions as illustrated in the left panel of Figure H] for 
the constituent non-strange quark mass. As a consequence, a first order phase transition is 
found at the critical chemical potential \f B = 361.7 MeV. The stable solutions are realized 
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FIG. 4. Left: grand potential for the stable, metastable and unstable phases, used to determine 
fji c £ and u quark mass; right: baryonic density as function of p,s- 

by the minimum of the thermodynamic potential. When stable and metastable solutions 
give the same value for the thermodynamic potential, the phase transition occurs as illus- 
trated in Figure 0] (left panel) . The phase of broken symmetry is realized for pb < p>% and 
the "symmetric" phase is realized for \xb > p%- At this crossing point of the curve, the 
two phases are in thermal and chemical equilibrium (Gibbs criteria). The baryon density, 
represented in Figure 0] (right panel) as function of ps, is given by the slope of the curve 
— Q as indicated by Equation (123]) . This quantity is plotted in Figured] (right panel) which 
shows that the condition of thermodynamic stability (d(—Q)/dpB > 0) is violated by the 
portion of the curve with negative curvature (unstable phase). 

This information can be complemented by the behavior of the pressure/energy per particle 
as a function of the baryonic density. To this purpose let us analyze the curve at T = 
in Figure [5] The pressure has three zeros that correspond to the extrema of the energy 
per particle. The third zero of the pressure, at ps = 2.36p 0) corresponds to an absolute 
minimum of the energy (see Figure [5] right panel). This is an important point of the model 
calculation, and the set of parameters is chosen in order to insure such a condition. In fact, 
the occurrence of an absolute minimum of the energy allows for the existence of finite droplets 
in mechanical equilibrium with the vacuum at zero-pressure (P = 0). For densities above 
a critical value, p c g = 2.36poj the system returns to a uniform gas phase. The equilibrium 
configuration for densities < ps < p% is> therefore, a mixed phase. 
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FIG. 5. Pressure (left) and energy per particle (right) as a function of Pb/po at different temper- 
atures (po = 0.17 fori -3 is the normal nuclear matter density). The points A and B (left panel) 
illustrate the Gibbs criteria. Only in the T = line the zero-pressure point is located at the 
minimum of the energy per particle. 



In view of the behavior above described, we can conclude that, for T = 0, the uniform 
non-zero density phase will break up into stable droplets, with zero pressure and density 
p% = 2.36po) in which chiral symmet ry is part ially restored, surrounded by a nontrivial 



vacuum with p# = P = (see also [36j, |107H111| ). In fact, for our choic e of the parameters 



the critical point at T = satisfies to the condition \x% < M^ ac 36j, Ill2[ |. where Mf ac is the 

mass of the z-quark in the vacuum. This can be seen by comparing pfg = 361.7 MeV (see 

the T-axis of Figure left panel) with the quark masses M v u ac = M v d ac = 367.7 MeV and 
M vac = 549 _ 5 MeV _ 



B. Phase Transition at Finite Temperature and Density/chemical Potential 

Let us now analyze what happens at T / 0. As the temperature increases (see Figure 
El left panel) the first order transition persists up to the CEP. At the CEP the chiral 
transition becomes a second order one. Along the line of a first order phase transition 
the thermodynamic potential has two degenerate minima, that are separated by a finite 
potential barrier making the potential non-convex. The height of the barrier decreases as 
the temperature increases and disappears at the CEP. Again, this pattern is characteristic 
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FIG. 6. Phase diagram in the SU(3) PNJL model. The left (right) part corresponds to the T — ps 
(T — pb) plane. Solid (dashed) line shows the location of the first order (crossover) transition. The 
dashed lines shows the location of the spinodal boundaries of the two phase transitions (shown by 
shading in the right plot). 

of a first order phase transition: the two minima correspond, respectively, to the phases 
of broken and restored symmetry. The borders of the coexistence area are marked by the 
dotted lines in Figure [6] which shows the T — fis phase diagram. The domain between the 
two dashed lines has metastable states which are characterized by large fluctuations. They 
are also solutions of the gap equations but their thermodynamic potential is higher than 
for the stable solutions. The left dashed curves represent the beginning of the metastable 
solutions of restored symmetry in the phase of broken symmetry, while the right dashed 
curves represent the end of the metastable solutions of broken symmetry in the restored 
symmetric phase. We also represent in Figure [6] (right panel) the region where the solutions 
of the gap equations are unstable. 

The location of the CEP is found to be at T CEP = 155.80 MeV and p c B EP = 1.87p 
{^ EP = 290.67 MeV). For temperatures above the CEP the thermodynamic potential has 
only one minimum and the transition is washed out: a smooth crossover takes place. In the 
left panel of Figure E] the crossover is represented by the dashed-dotted curve. To calculate 
the crossover we took the zero of the d 2 (q u q u )/dT 2 , i.e. the inflection point of the quark 
condensate (q u q u ). The crossover is defined as the rapid decrease of the quark condensate 
(the order-like parameter). To determine the range of the crossover region we choose to 
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define it as the interval between the two peaks around the zero of d 2 {q u q u )/dT 2 . This area 
is presented in gray in Figure [6] (left panel) . We can see that as //# increases the area where 
the crossover takes place is narrowed until we reaches the CEP. 

Finally, we will focus again on the energy per baryon. In Figure [5] (right panel), we plot 
the density dependence of the energy per baryon at different temperatures. We observe that 
the two points (zero of the pressure and minimum of the energy density) are not the same 
at finite temperature. In fact, as can be seen from Figure [5] (left panel), states with zero 
pressure are only possible up to the maximal temperature T m ~ 111 MeV. For T < T m the 
zero-pressure states are in the metastable density region and, as soon as T ^ 0, they do not 
coincide with the minimum of the energy per particle (see Figure [5] right panel). 

The arguments just presented allow to clarify the difference between confined quark 
matter (in hadrons) and bounded quark matter (droplets of quarks). 

This pattern of phase transition is similar to the liquid-gas transition in nuclear matter, 
a consequence of the fact that nuclear matter assumes its ground state at a non-vanishing 
baryon density ps ~ 0.17 fm~ 3 when T = 0. 

VI. NERNST PRINCIPLE AND ISENTROPIC TRAJECTORIES 



The isentropic lines contain important information about the conditions that are supposed 
to be realized in heavy ion collisions. Most of the stu dies on this topic have been done with 
lattice calculations for t wo flavor QCD at finite p [113| but there are also studies using 



different type of models 



112 



114 



1151 ] . Some model calculations predict that in a region 



around the CEP the properties of matter are only slowly modifie d as the collision energy is 
changed, as a consequence of the attractor character of the CEP 116 ]. 

Our numerical results for the isentropic lines in the (T, (Xb) plane are shown in Figured 
We start the discussion by analyzing the behavior of the isentropic lines in the limit T — > 0. 
As already analyzed in Section [V] our convenient choice of the model parameters allows a 
better description of the first order transition than other treatments of the NJL (PNJL) 
model. This choice is c rucial to obtain important results: the criterion of stability of the 
quark droplets 36|, Il09| is fulfilled, and, in addition, simple thermodynamic expectations 
in the limit T — > are verified. In fact, in this limit s — > according to the third law of 
thermodynamics and, as ps — > too, the satisfaction of the condition s/ p B = const, is 
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FIG. 7. Isentropic trajectories in the (T, ps) plane. The following values of the entropy per baryon 
number have been considered: s/ps = 1, 2, 3, 4, 5, 6, 8, 10, 15, 20 (anticlockwise direction). 

insured. We recall (Section [V} that, at T = 0, we are in the presence of droplets (states in 
mechanical equilibrium with the vacuum state (ps = 0) at P = 0). 

At T ytz 0, in the first order line , the behavior we find is somewhat different from those 
claimed by other authors 114J Ill7| ] where a phenomena of focusing of trajectories towards 
the CEP is observed. We see that the isentropic lines with s/ps = 1,...,6 come from the 
region of symmetry partially restored and attain directly the phase transition; the trajectory 
s /Pb — 1 goes along with the phase transition as T decreases until it reaches T = 0; and 
the other trajectories enter the hadronic phase where the symmetry is still broken and, after 
that, also converge to the horizontal axes (T = 0). Consequently, even without reheating in 



the mixed phase as verified in the "zigzag" shape of 113M115 



118 



all isentropic trajectories 



directly terminate in the end of the first order transition line at T = 0. 



In the crossover region the beha vior of the isentropic lin es is 



one obtained in lattice calculations 113| or in some models 114 



qua 



115 



itatively similar to the 



1191 ] . The trajectories 



with s/ Pb > 6 go directly to the crossover region and display a smooth behavior, although 
those that pass in the neighborhood of the CEP show a slightly kink behavior. 

In conclusion, all the trajectories directly term inate in the same point of the horizontal 
axes at T = 0. As already pointed out in 1121 ] . the picture provided here is a natural 
result in these type of quark models with no change in the number of degrees of freedom of 
the system in the two phases. As the temperature decreases a first order phase transition 
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occurs, the latent heat increases and the formation of the mixed phase is thermodynamically 
favored. 

We point out again that, in the limit T — > 0, it is verified t hat s — > and pb — > 0, as it 



should be. This behavior is in contrast to the one reported in 112| (see right panel of Fig. 
9 therein), where the NJL model in the SU(2) sector is used. The difference is due to our 
more convenient choice of the model parameters, mainly a lower value of the cutoff. This 
can be explained by the presence of droplets at T = whose stability is quite sensitive to 
the choice of the model parameters. 

VII. EFFECTS OF STRANGENESS AND ANOMALY STRENGTH ON THE 
CRITICAL END POINT 

As already noticed, an important question which emerges in the study of the QCD phase 
diagram is whether the hadronic phase and the quark gluon like phase are separated by a 
phase transition in the thermodynamical sense. The identification of possible critical points 
associated to the phase transition is a challenging problem. Since we are working within a 
model with (2+1) flavors with the 11^(1) anomaly included, we can explore new scenarios 
with regards to the usual SU(2) model without explicit anomaly. We have now two more 
parameters: the strange current quark mass, m s , and the anomaly coupling strength, g^, 
and it is certainly interesting to study the effect on thermodynamic quantities of changing 
these parameters. The anomaly is present via the 't Hooft interaction and its effects appear 
explicitly in the gap equations f fT8|) and in the expression of the thermodynamical poten- 
tial ( 1T9l through products of the anomaly coefficient by quark condensates. Here we will 
discuss the influence on the location and the nature of the critical points of the degree of 
explicit symmetry breaking of both chiral symmetry (in the strange sector) and axial Ua(1) 
symmetry. 

A. Role of the Strangeness in the Location of the CEP/TCP 

Let us study the influence of the strangeness in the location of the CEP/TCP. As we 
have already seen, the location of the CEP is found at T CEP = 155.80 MeV and ji % E1 



290.67 MeV (pb EP = l-87p ) when we use physical values of the quark masses [95|, Il20| : 
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FIG. 8. Left panel: the phase diagram in the SU(3) PNJL model. The solid lines represent the first 
order phase transition, the dotted line the second order phase transition, and the dashed line the 
crossover transition. Right panel: the phase diagram and the "line" of TCPs for different values 
of m s (the dotted lines are just drawn to guide the eye); the TCPs in both figures are obtained in 
the limit m u = rrid = and m s ^ 0. 



m u = m d = 5.5 MeV, m s = 140.7 MeV. 

A first point to be noticed is that in the PNJL model, contrarily to what happens in 
;he chiral limit only for the SU(2) sector (m u = m^ = 0, m s 7-- 0) where the TCP is found 
76|, when the total chiral limit is considered (m u = m^ = m s = 0), the phase diagram does 



not exhibit a TCP: chiral symmetry is restored via a first order transition for all baryonic 
chemical potentials and temperatures (see left panel of Figure [H])- Both situations are in 
agreement with what is expected: the chiral pha se transition at the chiral limit is of second 
order for Nf = 2 and first order for Nf > 3 121 ]. 

To study the influence of strangeness on the location of the critical points, we vary the 
current quark mass m s , keeping the SU(2) sector in the chiral limit and the other model 
parameters fixed. The phase diagram is presented in Figure M (right panel) as a function of 
fiB and T, and we consider different cases for the current quark mass m s . 

The pattern of chiral sym metry restoration via first order phase transition remains for 



m„ 



m d 



and m s < m 



crit 



122J. The va 



in lattice 123| or in model calculations 122 



or mf is a subject of debate; those found 
124| are lower than the physical strange current 
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quark mass (m s ~ 150 MeV). We found mf lt ~ 9 MeV in our model, low er th an lattice 



values 123l | and half of the value obtaine d in N JL model (mf lt = 18.3 MeV 



1201), but still 



consistent with other models of this type 124]. When m s > mf lt , at \xb = 0, the transition 



is of the second order and, as hb increases, the line of the second order phase transition will 
end in a first order line at the TCP. Several TCPs are plotted for different values of m s in 
the right panel of Figure As m s increases, the value of T for this "line" of TCPs decreases 
as hb increases getting closer to the CEP and, when m s = 140.7 MeV, it starts to move 
away from the CEP. The TCP for m s = 140.7 MeV is the closest to the CEP and is located 
at jj,l CP = 206.95 MeV and T TCP = 187.83 MeV. If we choose m u = m d ^ 0, instead of a 
second order transition we have a smooth crossover for all the values of m s and the "line" 
of TCPs becomes a "line" of CEPs. 

B. Role of the Anomaly Strength in the Location of the CEP 

The axial Ua(1) symmetry is broken explicitly by instantons, leaving a SU(N/)x SU(Nj) 
symmetry which determines the chiral dynamics. Since instantons are screened in a hot or 
dense environment, the Ua(1) symmetry may be effectively restored in matter. So, it is 
instructive to review this enlargement of chiral symmetry to SU(Nj)x SU(Nj)xUyi(l) in 
more detail. 

As already referred, the location and even the existence of the CEP in the phase diagram 



is a matter of debate 



51] . While different lattice calculations predict the existence of a CEP 



125 L the absence of the CEP in the phase diagram was seen in recent lattice QCD results 



79h81| , where the first order phase transition region near [ib = shrinks in the quark mass 



and fiB space when /j,b is increased |79l-l8l]. Due to the importance of the Ua(1) anomaly, 
already emphasized in Section [Til an d its influence on several observables, it is demanding 
to investigate possible changes in the location of the CEP in the (T, //#) plane when the 
anomaly strength is modified. In Figure [9] we show the location of the CEP for several 
values of go compared to the results for gr, , the value used for the vacuum. As already 
pointed out by K. Fukushima in [70], we also observe that the location of the CEP depends 
on the value of go and, when gr, is about 50% of its value in the vacuum, the QCD critical 
point disappears from the phase diagram; the first order region becomes narrower with the 
decreasing of the go value. 
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The results show that, in the framework of this model, the existence or not of the CEP 
is determined by the strength of the anomaly coupling, the CEP getting closer to the /x# 
axis as go decreases. 

VIII. SUSCEPTIBILITIES AND CRITICAL BEHAVIOR IN THE VICINITY OF 
THE CEP 



In the last years, the phenomenological relevance of fluctuations in the finite temperature 
and chemical p otent ial around the CEP/TCP of QCD has been attracting the attention of 
several authors 1261 ] . As a matter of fact, fluctuations are supposed to represent signatures 
of phase transitions of strongly interacting matter. In particular, the quark number suscep- 
tibility plays a role in the calculation of event-by-event fluctuations of conserved quantities 
such as the net baryon number. Across the quark hadron phase transition they are expected 
to become large; that can be interpreted as an indication for a critical behavior. We also 
remember the important role of the second derivative of the pressure for second order points 
like the CEP. 

The grand canonical potential (or the pressure) contains the relevant information on 
thermodynamic bulk properties of a medium. Susceptibilities, being second order derivatives 
of the pressure in both chemical potential and temperature directions, are related to those 
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FIG. 10. Phase diagram: the size of the critical region is plotted for Xb/Xb^ = 1(2)- The TCP 
is found for m u = rrid = MeV and m s = 140.7 MeV. 



fluctuations. The relevance of these physical observables is related with the size of the 
critical region around the CEP which can be found by calculating the specific heat, the 
baryon number susceptibility, and their critical behaviors. The si ze o f this critical region is 
important for future searches of the CEP in heavy ion-collisions 114 1 . 

The way to estimate the critical region around the CEP is to calculate the dimensionless 
ratio Xb/Xb ^ where Xb^ ^ s ^ e chiral susceptibility of a free massless quark gas. Figure ITOl 
shows a contour plot for two fixed ratios Xb/Xb^ = 1-0 anc ^ 2.0 in the phase diagram around 
the CEP. In the direction parallel to the first order transition line and to the crossover, it 
can be seen an elongation of the region where xb is enhanced, indicating that the critical 
region is heavily stretched in that direction. It means that the divergence of the correlation 
length at the CEP affects the phase diagram quite far fro m th e CEP and a careful analysis 
including effects beyond the mean-field needs to be done 127 ]. 

One of the main effects of the Polyakov loop is to shorten the temperature range where 



the crossover occurs 



44J. On the other hand, this behavior is boosted by the choice of the 



regularization (A — > oo) 76]. The combination of both effects results in high er b aryonic 
susceptibilities even far from the CEP when compared with the NJL model 128| . This 
effect of the Polyakov loop is driven by the fact that the one- and two-quark Boltzmann 
factors are controlled by a factor proportional to $: at small temperature $ ~ results 
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Specific heat as a function of T for different values of ^b around the CEP: [jf B rEP = 290.67 MeV 
and (i B = ^ C B EP ± 10 MeV. 

in a suppression of these contributions (see Equation (TT9|) ) leading to a partial restoration 
of the color symmetry. Indeed, the fact that only the 3— quark Boltzmann factors e 3 ^ Ep 
contribute to the thermodynamical potential at low temperature, may be interpreted as the 
production of a thermal bath containing only colorless 3-quark contributions. When the 
temperature increases, $ goes quickly to 1 resulting in a (partial) restoration of the chiral 
symmetry occurring in a shorter temperature range. The crossover taking place in a smaller 
T range can be interpreted as a crossover transition closest to a second order one. This 
"faster" crossover may explain the elongation of the critical region giving raise to a greater 
correlation length even far from the CEP. 

Now we show in Figure (TT] the behavior of the baryon n umber suscept ibility, xb, and the 



126 



129 ]. 



specific heat, C, at the CEP, which is in accordance with [120 . 

The baryon number susceptibility is plotted (left panel) as a function of the baryon 
chemical potential for three values of the temperature. The divergent behavior of xb at 
the CEP is an indication of the second order phase transition at this point. The curve for 
T > T CEP corresponds to the crossover and the other to the first order transition. 

Now we will pay attention to the specific heat (Equation ( 127)) ) which is plotted as a 
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function of the temperature around the CEP (right panel). The divergent behavior is the 
signal of the location of the CEP. The curve fj,B < ^ C b EP corresponds to the crossover and 
the other to the first order transition. 

The large enhancement of the baryon number susceptibility and the specific heat at the 
CEP may be used as a signal of the existence and identification of phase transitions in the 
quark matter. 

IX. CONCLUSIONS 

The concept of symmetries is a very important topic in physics. It has given a fruitful 
insight into the relationships between different areas, and has contributed to the unification 
of several phenomena, as shown in recent achievements in nuclear and hadronic physics. 

Critical phenomena in hot QCD have been studied in the framework of the PNJL model, 
as an important issue to determine the order of the chiral phase transition as a function of the 
temperature and the quark chemical potential. Symmetry arguments show that the phase 
transition should be a first order one in the chiral limit (m u = m^ = m s = 0). Working out 
of the chiral limit, at which both chiral and center symmetries are explicitly broken, a CEP 
which separates first and crossover lines is found, and the corresponding order parameters 
are analyzed. 

The sets of parameters used is compatible with the formation of stable droplets at zero 
temperature, insuring the satisfaction of important thermodynamic expectations like the 
Nernst principle. Other important role is played by the regularization procedure which, by 
allowing high momentum quark states, is essential to obtain the required increase of extensive 
thermodynamic quantities, insuring the convergence to the Stefan-Boltzmann (SB) limit of 
QCD. In this context the gluonic degrees of freedom also play a special role. 

We also discussed the effect of the Ua(1) axial symmetry both at zero and at finite 
temperature. We analyzed the effect of the anomalous coupling strength in the location of 
the CEP. We proved that the location of the CEP depends on the value of go and, when go 
is about 50% of its value in the vacuum, the QCD critical point disappears from the phase 
diagram. One expects that, above a certain critical temperature T e ff, the chiral and axial 
symmetries will be effectively restored. The behavior of some given observables signals the 
effective restoration of these symmetries: for instance, the topological susceptibility vanishes. 
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The successful comparison with lattice results shows that the model calculation provides 
a convenient tool to obtain information for systems from zero to non-zero chemical potential 
which is of particular importance for the knowledge of the equation of state of hot and 
dense matter. Although the results here presented relies on the chiral/deconfmement phase 
transition, the relevant physics involved is also useful to understand other phase transitions 
sharing similar features. 
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APPENDIX 

The quark condensates (<Mj), with i,j, k = u,d,s (to be fixed in cyclic order), determined 
in a self-consistent way, are given by: 

fe> = -2N c j^^-[9(A 2 -f)-fi + \E l )-f^\E l )} (29) 



where Ei is the quasi-particle energy for the quark i: Ei = y^p 2 + Mf. 

All calculations in the NJL model can be generalized to the PNJL one, as was shown 



in 44], by introducing the modified Fermi-Dirac distribution functions for particles and 



antiparticles used in the expression of (qiqi): 



<§) e -PEp J r 2<&e~ 2 P Ep -\-e~ 3 P Ep 



, v $ e -/3< + 2$e" 2/3 ^ + e~ 3l5E p 
ft\E t ) = * n (31) 

The mean field equations are obtained by minimizing the thermodynamical potential 
density ( 1T9J1 with respect to (<&<&) (i = u, d, s), $, and $. The respective mean field equations 
are the already given gap equation (fTBl) and 
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